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We consider changes of the topological charge of vortices 
in quantum mechanics by investigating analytical examples 
where the creation or annihilation of vortices occurs. In 
classical hydrodynamics of non-viscous fluids the Helmholtz- 
Kelvin theorem ensures that the velocity field circulation is 
conserved. We discuss applicability of the theorem in the 
hydrodynamical formulation of quantum mechanics showing 
that the assumptions of the theorem may be broken in quan- 
tum evolution of the wavefunction leading to a change of the 
topological charge. 



I. INTRODUCTION 

Vortices can be found both in classical and quantum 
physics. One can encounter vortices, for instance, in wa- 
ter that spins around or in the air as a ring of smoke 
[1]. Quantum mechanics can be formulated in a language 
of hydrodynamics (see e.g. [2,3] and references therein). 
Such a formulation (very useful also in quantum chem- 
istry [2]) provides a basis for definition of topological de- 
fects like vortices [4] , whose features are even more strik- 
ing than those that we find in classical physics. Quantum 
vortices appear not only in systems described by the lin- 
ear Schrodinger equation. Indeed, they were experimen- 
tally observed in superfluid Hell [5] and a Bose-Einstein 
condensate (BEC) of trapped alkali atoms [6,7] which 
are typically described, within the mean-field approxi- 
mation, by a nonlinear equation [8]. 

In hydrodynamics of ordinary non-viscous fluids, the 
circulation of the velocity field is conserved in time evo- 
lution due to the celebrated Helmholtz-Kelvin theorem 
(HKT) [9]. This theorem is often employed in quantum 
mechanics [3,4,10,11]. However, uncritical usage of the 
HKT, may lead to the incorrect conclusion that stabil- 
ity of vortices and vortex rings in quantum fluids is fully 
guaranteed by the HKT [11]. In the present paper we 
discuss the basic assumptions of the HKT and show an- 
alytical examples where these assumptions can be easily 
broken by quantum evolution of the wave function, lead- 
ing to changes of the topological charge. Difficulties in 
fulfilling the assumptions of the HKT have been already 
pointed out in conclusions of Ref. [12]. 



II. HYDRODYNAMICAL FORMULATION OF 
QUANTUM MECHANICS 

To establish connections between quantum mechanics 
and fluid dynamics we write the wave function in the form 
\l/(f, t) = \J p(f, t) exp(ix(r, t)) [2], (where p(r,t) stands 
for density of a probability fluid) and define the velocity 
field 



-V X (r,i). 



(1) 



Provided v and its partial derivatives are well defined, we 
may rewrite the Schrodinger equation [with a potential 
V(r, t)] in the form 



dv 
m— + V 
at 



-mv + V 



h 2 v 2 Vp 

2m y/p 



= 0, 



dp 
dt 



+ V ■ (vp) = 0. 



(2) 



(3) 



Equation (3) is an ordinary continuity equation, while 
Eq. (2), in the limit of H — > 0, becomes similar to the 
dynamical equation of a curl-free non- viscous fluid [9] . 

In quantum mechanics the wave function has to be 
single valued. To satisfy this condition one arrives at 
the Feynman-Onsager quantization condition [13] for the 
circulation of the velocity field Tc around any closed con- 
tour C 



v ■ dl = 



(4) 



where n — 0, ±1,±2, . . .. Due to the definition (1), we 
may not expect n ^ unless at a certain point on a sur- 
face encircled by the contour C, V x v becomes singular 
(behaves as a Dirac delta function, see for example [5]). 
We refer to n as a topological charge since any continuous 
deformation of the contour C, which does not incorpo- 
rate any other points where curl of the velocity field does 
not vanish, can not change Tc- 



III. THE HELMHOLTZ-KELVIN THEOREM 

One may ask whether the circulation Tc is conserved 
in time evolution of the velocity field. For a fixed contour, 
a possibly moving vortex may leave an area bounded by 
the contour and consequently Tc changes its value. Thus, 
the relevant situation we should consider corresponds to 



the case where we let an initially defined contour evolve 
in time according to the velocity field. 

Parameterizing the initial contour by £ variable, i.e. 
C(*o) = to)}, and employing the equation of motion 
of a contour 



*(a;,y) oc (x + iay)e-^ 2+Xy ^' 2 , 



(9) 



d_ 
~dt 



yields 



c(t) {^ +{ "-^ + ¥ v2 



(5) 



• dl. (6) 



If the contour is drawn through points where the velocity 
field v and its partial derivatives are well defined, apply- 
ing: Eq. (2), equality (v ■ V)v = \Vv 2 - v x (V x v), and 

using the fact that, at least on the contour V x v = 0, 
we get 



dt 1 



C(t) 



0. 



(7) 



Equation (7) establishes the Hchnholtz-Kelvin theorem 
originally proved for a non- viscous fluid [9]. Applicabil- 
ity of this theorem to a non-linear quantum systems is 
presented in the Appendix. 

The HKT is fulfilled provided the contour C(t) that 
goes initially through points where the velocity field is 
well defined will evolve in time through such points only. 
In the following we will see how such an assumption may 
be broken by quantum evolution of the wave function 
leading to a change of the topological charge. 



IV. ANALYTICAL EXAMPLES 

In this section we discuss a few possibilities of breaking 
of the HKT. To this end we analyze three different an- 
alytically solvable examples where breaking of the HKT 
becomes evident. 

First example shows that changes of the topological 
charge of a vortex correspond to the appearance of a 
nodal line of the wave function. In Ref. [14] the authors 
show the instability of a vortex placed in an anisotropic 
two-dimensional harmonic trap. We would like to discuss 
the source of violation of the HKT in such a case. We 
will calculate the velocity field and analyze changes of its 
circulation from a point of view of the HKT. 

Consider an anisotropic two-dimensional harmonic os- 
cillator with the potential 



1 A 2 

V(x,y) = -x 2 + —y\ 



(8) 



(where we use the units of the harmonic oscillator corre- 
sponding to the x-direction). The initial wave function 
is prepared as a superposition of the two lowest excited 
states 



with a real parameter a > 0. Choosing any contour 
C that encircles the origin of the coordinate frame we 
find out that the circulation Tc corresponds to n = 1 
vortex. However, the circulation around that point is 
not conserved in time. Indeed, analysis of time evolution 
of the velocity field 



v(x,y,t) 



a cos(Et) 



x 2 + y 2 a 2 + 2axy sin(Et) 



(~ye x + xe y 



(10) 



(where E = A — 1 is the energy difference of the two 
lowest excited states) reveals that for Et <G [0, 7r/2), Tc 
corresponds to n = 1 while for Et € (tt/2, it] to n = —1. 

To discuss applicability of the HKT we have to in- 
vestigate time evolution of the initially defined contour. 
The problem requires careful treatment because for Et = 
tt/2 — e (where e — > + ) the velocity field goes to zero 
everywhere except on the nodal line, x + ay = 0, where 
it diverges v — > j^{e x + ae y ). To demonstrate violation 
of the HKT we have to show that a contour initially en- 
circling the center is not extended to infinity but indeed 
faces singularity for Et = w/2. 




FIG. 1. Plot of the contour encircling the center of the trap 
for Et — » tt/2 (dashed line). The contour was initially (t — 0) 
chosen to be a circle with the unit radius centered at (0.5, 0) 
(solid line). The thick solid line represents the nodal line 
given by the equation x + ay — 0. Initially, the point of the 
contour closest to the origin of the coordinates is situated at 
the distance 0.5. Therefore in course of the evolution none of 
the points of the contour can get closer to the origin than at 
the distance 0.5, i.e. can not cross the dotted circle with the 
radius 0.5 that is centered at (0,0). For this presentation we 
have chosen A = Vz, a = 1 — see text for details. 

Suppose we define a contour, encircling the origin, at 
a time t = 0, and look for a time evolution of its arbi- 



trary point i, whose position we denote by {xi(t), yi(t)}. 
Combining Eq. (5) and Eq. (10) we get 

^ = -f{xi,Vi) yu ^ = f(xi,yi) Xi, (11) 



with 



f(x, y) = 



acos(Et) 



x 2 + y 2 a 2 + 2axy sm(Et) 



(12) 



From Eq. (11) it is apparent that the distance of the point 
{xi(t), yi(t)} from the origin does not change in course 
of the evolution. Indeed the quantity Xi(t) 2 + yi(t) 2 is a 
constant of motion as long as the function (12) is well de- 
fined. Note that it diverges, at Et — ir/2, for points lying 
on the nodal line. Therefore, for Et = ir/2 — e with ar- 
bitrary small e > nothing extraordinary happens with 
the contour until at Et = it/ 2 the contour faces singular- 
ity (see Fig. 1 for an example of the contour evolution). 
Thus, the HKT can not be applied because quantum evo- 
lution of the wave function pushes an evolving contour 
to singular points and the quantities involved in the the- 
orem become undefined. 

We would like to stress that the presented scenario 
where the vortex situated in the harmonic trap reveals 
periodic changes of the topological charge is possible only 
when the initial vortex is placed exactly at the center 
of the potential. Otherwise the vortex changes its po- 
sition — it escapes to infinity and another vortex with 
the opposite charge arrives from infinity and the scenario 
repeats periodically. 

In the second example we would like to present vio- 
lation of the HKT due to annihilation of a vortex ring. 
This example comes from Ref. [4] where creation and 
annihilation of a vortex ring for freely moving particle is 
presented. Even though it is not clearly stated in Ref. [4], 
the HKT is not fulfilled in that case. The wave function 
of a vortex ring [4] (in dimcnsionless units) may be writ- 
ten in the form 



i&(x,y,z,t) oc [(x - kt) 2 



z 2 -l 



+ ^3(z + ^)]e ^fe^; - ^fc2t / 2 , 



(13) 



where the wave vector k = (k, 0, 0) is related to a motion 
of the "center of mass" of the vortex ring with a constant 
velocity k. The vortex ring corresponds to the nodal line 
of the wave function and is located at the intersection of 
the plane z + t = and sphere (x — kt) 2 + y 2 + z 2 = 1. 
At time t = —1, the vortex is born at a point (— k, 0, 1) 
and, at time t = 1, it disappears at a point (k, 0,-1). 
The radius of the vortex changes in time as \f\ — t 2 . 

Suppose at any time t € (— 1, 1), we define a contour 
C so that it encircles the vortex (see Fig. 2) and the 
circulation of the velocity field corresponds ton = 1 . The 
evolving contour can not cross the vortex ring without 
facing a singularity in the velocity field. However, the 
ring at some moment starts shrinking and at t = 1 it 



reduces to a point. Consequently, at t = 1, the contour 
must go through a singularity of the velocity field and 
the integral (6) and also the HKT become meaningless. 




FIG. 2. Schematic plot of the nodal line corresponding to 
a vortex ring (solid line) and a contour that encircles the ring 
(dashed line). In last stage of time evolution the ring shrinks 
to a point and the contour faces singularity. 

The last example we would like to comment on shows 
collision of two vortices leading to the appearance of a 
single doubly charged vortex. Consider, for simplicity, 
a two-dimensional H atom initially in the first excited 
state with angular momentum L z = 1 [i.e. rpi(r,ip) oc 
rexp(— 2r/3) exp(ico) in atomic units] that is driven res- 
onantly by a circularly polarized electromagnetic field. 
The field frequency is tuned to the transition between 
the first and second excited energy eigenstates. It results 
in the familiar Rabi oscillation [15] between the ipi(r, ip) 
state and the second excited state with angular momen- 
tum L z — 2 [i.e. i/)2{r,ip) oc r 2 exp(— 2r/5) exp(i2<p)]. 
Indeed, the time evolution of the wavefunction reads 

|V(t)) oc cos(Dt)e- lElt + e- lE2t sm(Dt) \ip 2 ) , (14) 

where D is a dipole matrix element [15]; E\, Ei are en- 
ergies of l^i), IV^} states. 

For t = there is a vortex with n = 1 at the center 
of the coordinates as expected for the state with angular 
momentum L z = 1. When Dt increases another vortex 
with n = 1 moves in from infinity, collides with the first 
one (situated at the center of the coordinates during the 
whole time evolution) at Dt — ir/2 and then moves out 
to infinity again and so on. During the collision, i.e. for 
Dt = tt/2, a single vortex with n — 2 is formed. In the 
considered example, the HKT does not apply: if we de- 
fine a contour so that it encircles the vortex with n = 1 
situated at the center of the coordinate, such a contour 
encounters a singularity in the velocity field during the 
collision with the other vortex. Such a behavior, is ex- 
pected to hold in an arbitrary collision between two vor- 
tices leading to formation of one, doubly charged, vortex. 

These examples are not the only ones illustrating that 
vortices can disappear or change their charge in the 
course of time evolution. They were chosen to illustrate 
that such processes can happen in different physical sys- 
tems and results in the violation of the HKT. The theo- 
rem is explored in the context of vortices' stability (see 
e.g. [11]) to imply the constancy of a vortex topological 
charge. As we have clearly illustrated the HKT can not 
assure persistence of vortex currents in the quantum me- 
chanical systems. Therefore we conclude that the vortex 



topological charge is not as robust quantity as it is com- 
monly believed (see e. g. [14] and references therein, [16]). 
Finally, we would like to mention that appearance of vor- 
tices can happen in a way that does not violate the HKT. 
Namely, they can appear in the form of a closed vortex 
line that springs from a point or as a vortex-antivortex 
pair creation from a node of a wave function [4] . 



V. SUMMARY AND CONCLUSIONS 

We have analyzed the applicability of the Helmholtz- 
Kclvin theorem in the hydrodynamical formulation of the 
quantum mechanics. The velocity field of the probability 
fluid is defined as the gradient of the phase of the quan- 
tum wave function. This implies that nonzero circula- 
tion, along a given contour, may come out in the system 
if the field reveals a singularity at certain points on a 
surface encircled by the contour. Adopting the HKT to 
the quantum liquid may suggest that such a topological 
charge of the system can not change. 

However, the HKT may be employed if a given contour 
evolves through points where the velocity field is well de- 
fined. It may happen even in classical hydrodynamics 
that such an assumption is not fulfilled if, e.g., liquid 
encounters obstacles in the flow [9]. In quantum liquids 
the situation is more complicated because a singularity 
is necessary for a nonzero circulation. Indeed, phase of a 
wave function is undefined at a vortex core, which means 
that singularities appear in hydrodynamical formulation 
of quantum mechanics whenever vortices show up. This 
property makes distinction of hydrodynamical descrip- 
tion of quantum system from the classical hydrodynam- 
ics. We have presented simple analytical examples where 
the quantum evolution of a wave function pushes a con- 
tour to a singular point. Such a process is accompanied 
by a change of the vortex topological charge. 

We have illustrated the violation of the HKT choos- 
ing examples that correspond to linear quantum systems. 
Macroscopic quantum behavior is present in interact- 
ing many particle systems that are usually described in 
the mean field approximation by a nonlinear Schrodingcr 
equation. In Appendix we complete the analysis of the 
HKT for a nonlinear Schrodingcr evolution. Vortices 
maybe also investigated in propagation of classical light 
[16,17]. Indeed, the evolution equation of the slowly vary- 
ing envelope of a light beam can have identical form as 
the Schrodinger equation, which allows for exactly the 
same considerations as we have presented above. The 
dynamics of vortices and escapes of the off-center vortex 
in an analog of the asymmetric harmonic potential have 
been experimentally observed in such systems [18]. 
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cal reading of the manuscript. Support of KBN under 
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VI. APPENDIX: HELMHOLTZ-KELVIN 
THEOREM FOR NONLINEAR SYSTEMS 

For completeness we analyze here the application of the 
HKT in hydrodynamical formulation of nonlinear quan- 
tum mechanics. The so-called Gross-Pitaevskii equa- 
tion is successfully used to describe the properties of a 
Bose-Einstein condensate in trapped alkali atoms [8] . In 
particular this equation exhibits various vortex solutions 
[11]. To make the discussion more general we will also 
include the mean field term describing a possible dipolar 
interactions between condensed atoms [19]. The nonlin- 
ear Schrodinger equation of interest reads 

d fi 2 
ih—V{r,t) = - —V 2 y(r,t) + V(r,t)V(r,t) 



+ 



J d 3 r'V mt (r-7)\*(7,t)\ 2 ^(r,t), (15) 



where Vi nt corresponds to two body interactions (e.g. 
point or dipolar interactions [8,19]). Substitution 
^(r, t) = p(r, t) exp(ix(r 7 1)) leads to the continuity 
equation (3) and the following dynamical equation for 
the velocity field 



rn 



dv 

di ' 



vQmu 2 + V + J d 3 r'V lnt (r-f')p(?,t) 



n 2 v 2 yp 

2m sj~p 



= 0. 



(16) 



The only difference in comparison with Eq. (2) is the ap- 
pearance of two new terms in the bracket of Eq. (16). 
These terms, however, do not change the proof of the 
HKT — for example Eq. (6) does not possess any poten- 
tial like terms. Therefore the Helmholtz-Kelvin theorem 
holds also for systems described by Eq. (15) under the 
same assumptions as in the case of linear quantum me- 
chanics discussed in Sec. III. 
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